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We investigate the spin behavior of close-in rocky planets and the implications 
for their orbital evolution. Considering that the planet rotation evolves under simul- 
taneous actions of the torque due to the equatorial deformation and the tidal torque, 
both raised by the central star, we analyze the possibility of temporary captures in 
spin-orbit resonances. 

The results of the numerical simulations of the exact equations of motions indicate 
that, whenever the planet rotation is trapped in a resonant motion, the orbital decay 
and the eccentricity damping are faster than the ones in which the rotation follows 
the so-called pseudo-synchronization. Analytical results obtained through the aver- 
aged equations of the spin-orbit problem show a good agreement with the numerical 
simulations. 

We apply the analysis to the cases of the recently discovered hot super-Earths 
Kepler-lOb, GJ 3634b and 55 Cnce. The simulated dynamical history of these sys- 
tems indicates the possibility of capture in several spin-orbit resonances; particularly, 
GJ 3634b and 55 Cnce can currently evolve under a non-synchronous resonant mo- 
tion for suitable values of the parameters. Moreover, 55 Cnc e may avoid a chaotic 
rotation behavior by evolving towards synchronization through successive temporary 
resonant trappings. 

Key words: celestial mechanics - planets and satellites: general. 
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1 INTRODUCTION 

Tidal theories have gained great attention in the last decade due to the large number of discovered close-in planets which 
orbit their parent stars with periods of a few days. According to classical theories (e.g. Darwin 1880; Kaula 1964; Mignard 
1979), the long-term tidal evolution of these planets results in orbital decay and circularization, on time-scales which depend 
on the physical properties of the interacting systems (Hut 1981; Dobbs-Dixon, Lin & Mardling 2004; Ferraz-Mello, Rodriguez 
& Hussmann 2008; Jackson, Barnes & Greenberg 2009; Rodriguez & Ferraz-Mello 2010). 

The rotation of short-period planets is also affected by tidal interactions. It is nowadays accepted that the spin rate 
of planets soon after their formation was much higher than their orbital mean motion (e.g. Agnor, Canup & Levison 1999; 
Kokubo & Ida 2007). However, the cumulative tidal effects causes the planet to spin down, reaching its stationary value in 
time-scales much shorter than those needed to circularize the planet orbit. This evolution ultimately results in the synchronous 
rotation (the spin state in which the orbital and rotation periods coincide), which occurs when the planet orbit is circularized 
(Hut 1981; Ferraz-Mello et al. 2008). For non-circular orbits, the synchronous rotation only occurs if an additional torque acts 
to counteract the tidal torque. 

The additional gravitational torque may be generated by planets with a rocky-like composition having a permanent 
equatorial deformation or frozen-in figure (see Goldreich & Peale 1966; Greenberg & Weidenschilling 1984). These torques, in 
contrast to the tidal ones, are conservative. In the context of the two-body problem, the gravitational interaction of the star 
with an asymmetric planet drive the planet rotation into different regimes of motion, including oscillations around exact spin- 
orbit resonances (Goldreich & Peale 1966). When dissipative tidal effects are taken into account, the oscillations are damped 
and the planet rotation is trapped into a resonance. Indeed, the rotational motion of some Solar System bodies can be 
explained by the combined effects of the tidal and asymmetric torques (Goldreich & Peale 1966; Greenberg & Weidenschilling 
1984; Correia & Laskar 2009). 

Mercury is an example of a rocky planet with a rotation state currently evolving in the 3/2 spin-orbit resonance (Pettengill 
& Dyce 1965). Several models of the secular spin evolution of Mercury have been developed, including i) the chaotic evolution 
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of its orbit (Correia & Laskar 2004); ii) effects of the core- mantle friction (Correia 2006; Correia & Laskar 2009); iii) the 
evolution starting with retrograde and subsequent synchronous motion combined with impact events (Correia & Laskar, 2012; 
Wieczorek et al. 2012); iv) the consideration of non- linear tidal theories (Makarov 2012). However, despite the large amount 
of previous works, the question on how the current resonant state of Mercury was achieved still remains under discussion. 

Concerning the rotation of Earth-like close-in extrasolar planets, a near synchronous motion is usually assumed in view 
of tidal evolution. However, Correia, Levrard & Laskar (2008) have shown that when an additional torque is considered, the 
equilibrium rotation of tidally evolved planets can assume rotational states which depart from synchronism. Moreover, as the 
capture probability in a non- synchronous spin-orbit resonance requires a non-circular orbit (Goldreich & Peale 1966; Celletti 
& Chierchia 2008), the rotation of several recently discovered super-Earth planets in eccentric orbits can be currently evolving 
under a resonant motion, provided the planet figure can support a permanent equatorial deformation. As the eccentricity is 
tidally damped, several resonant trappings would become unstable and the planet rotation ultimately reaches synchronization 
for almost circular orbits. However, the locking in a non-synchronous resonance with the very short orbital period can protect 
certain planet's longitudes from solar radiation, providing more favorable conditions with respect to habitability as compared 
to non-resonant rotation (Dobrovolskis 2007). Moreover, Heller, Leconte & Barnes (2011) studied the equilibrium rotation of 
super-Earths and the implications on their habitability. They found that, under assumption of an evolution driven by tidal 
torques, the equilibrium value of the rotation period of GJ 581 d is probably about two times smaller than the orbital period. 

In a recent work, Callegari & Rodriguez (2012) have analysed the rotation of extrasolar planets with orbital periods less 
than 33 days and masses below the mass of Neptune. Depending on one parameter related to the equatorial ellipticity and 
also on the orbital eccentricity, they identified regions of chaotic rotation in the phase space of some particular planets like 
55 Cnce, GJ 1214 b, Gliese 876 d, among others. 

The present investigation aims at two main aspects: i) How the orbital evolution of the planet, due to tidal interaction 
with the central star, is altered when the asymmetric torque is considered; and ii) To analyse the occurrence of captures in 
spin-orbit resonant motions under combined effects of both torques for the planets Kepler-10 b, 55 Cnc e and GJ 3634 b. These 
planets are examples of hot super-Earths, which are close-in planets having a few Earth-masses with a rocky-like composition. 
We are interested in the past behavior in the history of the rotation and the subsequent evolution to synchronization. We 
also discuss the possibility of a current evolution of the rotation of these planets into a super-synchronous resonance state. 
Moreover, we investigate the implications on the orbital evolution due to such temporary resonant trappings. 

Despite the continuously growing number of discovered planets, their rotation state and equilibrium figures remain 
unknown. In this work, we perform several numerical simulations of the evolution of the systems, assuming several values of 
unknown parameters, such as equatorial ellipticity and initial rotation. The results allow us to understand the past history of 
these systems, their current orbital configuration and final rotational states. It should be noted that the majority of previous 
works concerning spin-orbit interactions with tidal dissipation have been addressed in the context of a fixed Keplerian orbit, 
with constant values of semi-major axis and eccentricity (e.g., Celletti, Froeschle & Lega 2007; Celletti & Chierchia 2008). 
This assumption is adequate for Solar System bodies, where tidal interactions are weak and, consequently, variations of orbital 
elements are small. 

On the other hand, the existence of extrasolar planets with orbital periods of a few days indicates that their current 
orbital configurations would be the consequence of strong tidal interaction with the host star. In this case, the Keplerian 
approximation is no longer valid. Thus, we analyse spin-orbit interactions as the orbit evolves due to tidal effects. Solving 
the exact equations of motion of the problem, we show that the forces arising from the potential of a asymmetric planet also 
contributes to the orbital variation, accelerating the orbital decay and the eccentricity damping, but only in the cases in which 
the rotation is locked in resonance. 

The spin-orbit problem can also be modeled by the use of averaged equations of motion, where only terms with secular 
and resonant arguments are retained (e.g. Celletti & Chierchia 2008; Correia & Laskar 2009). Through this model, we compute 
the averaged variations of semi-major axis and eccentricity due to the combination of asymmetric and tidal forces. We also 
give an explicit expression for the critical value of the equatorial prolateness, for which a given spin-orbit resonance would 
become unstable. The study of the particular cases of 2/1, 3/2 and 1/1 spin-orbit resonances indicates a good agreement with 
the results of the numerical simulations of the exact equations. 

This paper is structured as follows. In Section O we present the exact equations of motion of the spin-orbit problem. 
Section [3] shows the results of the numerical simulations, illustrating the evolutions of Kepler-10 b, GJ 3634 b and Cnc 55 e. 
We show that, assuming a nonzero value of the equatorial ellipticity, different temporary resonant trappings of the rotation 
can occur, with implications in the orbital decay and circularization. An analytical approach in the context of the averaged 
equations is shown in Section |4j explaining the results obtained in the numerical simulations. Finally, Section [5] is devoted to 
discussion and conclusions. 



2 THE MODEL 

We consider a planet orbiting a central star with a period of a few days. Hereafter, we refer to m and M as the planet and 
star masses, respectively; R and r as the planet radius and the instantaneous star-planet distance, respectively. The elements 
a,e,n stand for the semi-major axis, eccentricity and mean-motion of the astrocentric orbit of the planet, respectively. 

We suppose that the rotating planet body can be approximated by a homogeneous triaxial ellipsoid, with the equatorial 
axes a and b and the rotation axis c. We assume that a > b > c and c is normal to the orbital plane (i.e., no obliquity is 
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Figure 1. Schematic view of the ellipsoid, showing the angles used in the formulation of the problem. The long axis makes an angle 9 
with the reference axis X, which is assumed coincident with the orbital pericenter direction . The angle between the long axis and the 
star-planet direction is given by / — 8, where / is the true anomaly of the planet. 



considered; planar problem). In addition, we assume that the equatorial asymmetry (a 7^ b) corresponds to a frozen-in body or 
permanent deformation of the planet. The assumption above can be supported under consideration of a rocky composition of 
the planet body. This is the case of super-Earth planets, which are the subject of the present investigation. The gravitational 
potential of the planet at the distance r from its centre is given by (e.g. Beutler 2005, p. 102) 

U22 = 3^^C 2 2Cos2(/-0), (1) 

where G is the gravitational constant and the ang / and 6 are denned in Figure Q] Note that 6 = fi is the angular rotation 
velocity of the planet. The parameter C22 describes the equatorial ellipticity of the planet's gravity field and is related to the 
principal momenta of inertia A, B,C through C22 = (£/4) (B — A)/C, where £ is the structure constant given by £ = C / (rnR 2 ) 
and C is the polar moment of inertia. 

We neglect the contribution of terms arising from the oblateness of the rotating body, that is, J2 and higher order terms. 
In fact, in the planar problem, the gravitational potential corresponding to the J2-term depends only on the distance r. In 
this case, the potential generates a radial force with zero torque that does not affect the rotation evolution. 

At the distance r from the centre of the planet, the gravitational field U22 gives rise to a force on a point mass M which 
is given b}0 

.,9(722 , . M dU 2 2 » ,„s 

F22 = M VU22 = M— — rH — ; f , (2) 

or r of 

where f and f are unit vectors with origin in the star (see Fig. [!}. The force F22 alters the orbital component of the angular 
momentum, whereas the reaction force governs the evolution of the spin through the creation of the torque on the planet 
given by T22 = — r x F22 = —M(dU22/df) k, where r = r f and k = f x f . Thus, 

T22 = 6 Gm ^ R2 C22 sin 2(/ - 0)k. (3) 

Note that, for CCl = T22 and C22 = (£/4)(B — A)/C, the above equation agrees with the classical expression for the torque 
(e.g. Danby 1962; Goldreich & Peale 1966). 

The effects of the dynamical tides are raised by the gravitational action of the star on the inelastically deformed rotating 
planet. Due to the internal viscosity of the body, there exists a delay At (time lag) between the star perturbation and the 
maximum deformation. The relationship between At and tidal frequencies (which are linear combinations of and n) depends 
on the internal rheology, which is poorly known even for Solar System bodies. For sake of simplicity, in this work, we adopt 
a linear model with constant At (Darwin 1879; Mignard 1979; Hut 1981; Rodriguez et al. 2011a); for a review of other tidal 
models, the reader is referred to Efroimsky & Williams (2009) and Ferraz-Mello (2012). 

According to Mignard (1979), the tidal force and torque are given by 

F tidc = -3faAt - 10 [2r(r • v) + r 2 (r x + v)], (4) 



1 We choose the pericenter of the planet orbit as the origin for the calculation of the longitudes because the main contribution to the 
variation of the longitude of the pericenter comes from the zonal term proportional to J2, which has not been included in the present 
paper (see Ferraz-Mello et al. 2008). 

2 Since we use the expression for U22 adopted in Beutler (2005), the corresponding force is computed as the positive gradient of the 
potential. 
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Table 1. Physical parameters and current orbital elements of the analysed systems ("Batalha et al. 2011; b Gillon et al. 2012; c Bonfils 
et al. 2011). The first two planets were discovered through the transit technique, allowing us a good estimation of their radii, and 
Rgj stand for mass and radius of the Earth, respectively. 



System 


A/(m ) 


m(m ffi ) 




(^current ( a u) 


^current 


Kepler-10 a 


0.895 


4.56 


1.416 


0.01684 





55 Cnc 6 


0.905 


8.63 


2.17 


0.0156 


0.057 


GJ 3634 c 


0.45 


7.0 




0.0287 


0.08 



and 

T tidc = 3k 2 At GM r l R5 [-r 2 Cl + r x v], (5) 

where v = r and k 2 is the second degree Love number. In this paper, we only consider the tides raised by the star on the 
planet, neglecting the tides raised on the star by the planet. Indeed, the strength of stellar tide is proportional to the planet 
mass and can be safely neglected in the case of super-Earths (see Rodriguez & Ferraz-Mello 2010). Note however that, after 
circularization of the planet orbit, stellar tides may be the only source of orbital decay. 

In the reference frame centered in the star (see Figure [T]), the exact equations of motion of the planet can be now written 

as 

.. _ G(M + to) (M ■ 



Mm 



-(Ftide + Fa 



§ = n= ^(Ttide+Taa). (6) 

The perturbing forces F22 and F t id c alter the Keplerian orbit of the planet which is defined by the first term of the first 
equation, while the corresponding torques in the second equation produce variations of the planet rotation. Note that, in the 
planar case, this last equation can be written in scalar form, because both torques are normal to the orbital plane. 



3 NUMERICAL SIMULATIONS 

In this section we apply the analytical model described in the previous section, in order to investigate the motion of three 
known super-Earths planets, namely Kepler-10 b, 55 Cnc e and GJ 3634 b. For this task, the exact equations of motion ((6)l were 
integrated numerically using the RA15 code (Everhart 1985). We stress that no scaling parameter, often used to accelerate 
the tidal evolution, was used in our simulations (see, for instance, Tittemore & Wisdom (1988), for the validity of the usage 
of the scaling parameter). 

Table 1 shows the physical parameters and current orbital elements of the systems considered in this paper. 



3.1 Setting initial configurations 

Numerical integrations require a setting of the initial configurations of the systems. In the case of problems including dissipative 
processes, the choice is complicated, mainly, due to the inadequacy of the simulations of the past evolution of dissipative 
systems through integrations backward in time (see details in Michtchenko & Rodriguez 2011). However, as shown in that 
paper, the final state of a migrating system is independent of its initial position in the phase space, defined by three values 
of parameters: the mass and semi-major axes ratios and the total angular momentum. As a consequence, to assess the orbit 
configuration in the past, we have to set solely the starting semi-major axis value. The initial eccentricity should satisfy the 

conservation of the total angular momentum of the system. 

The total angular momentum of the system is composed of the orbital component written as to 
and the rotational component given by CO, = £mR 2 Q.. Using Kepler's third law, it is easy to show that the ratio of rotational 
and orbital components is of the order of (R/a) 2 (see Rodriguez et al. 2011a); since R <C a, the contribution of rotation to 
the total angular momentum can be neglected. Therefore, considering that the orbital angular momentum is nearly conserved 
during tidal migration, we obtain that 

eini — 4/I - (1 — e currcnt)j (?) 

where the subscripts "ini" and "current" stand for initial and current orbital elements, respectively. In this way, choosing an 
initial value for the semi-major axis and knowing the current values of a and e, equation (J7| allows us to constrain the initial 
value of the eccentricity (see next section and specially Figure [6}. It should be emphasized that the value obtained in such a 
way must be considered as an approximation of the initial value of the eccentricity, since the planets may belong to systems 
with two or more planets. 

The value of At is determined by the relationship between the tidal frequency v and the quality factor Q. In the case 
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Figure 2. The variation of the rotation with eccentricity of the Kepler-lOb orbit. The arrow indicates the direction of the eccentricity and 
rotation evolutions. The curve with segments line shows the stationary tidal solution (STS) of the rotation, indicating a good agreement 
with the result of numerical simulation for C22 = (solid black curve). The dotted line corresponds to the case C22 = 10 -5 , where 
the trappings 2/1 and 3/2 occur. After escaping from the 3/2 resonant motion, the rotation follows the STS towards a synchronization 
between the rotation and orbital motions. 



of a linear approximation, we have 1/Q = vAt, where At = const. It is known that, in the Fourier spectrum of the tidal 
potential, the most important frequencies (that is, those involved in terms of order zero and order one in the eccentricity) are 
^0 = 2f2 — 2n, v\ = 2fi — 3n and ^2 = 2S1 — n . The first one is associated to the term which is independent of the eccentricity, 
while the two last frequencies correspond to terms which are proportional to e (Darwin, 1879; Ferraz-Mello et al. 2008). Using 
the frequency vq and a modified version of the quality factor written in terms of the Love number &2 as Q' = 3Q/2/C2, we 
obtain 

k ^ = WW^y (8) 

Actually, the spin rotation of exoplanets is not available from observations, thus we choose its initial value arbitrarily but 
considering a super-synchronous motion (Q > n). 

Finally, according to Callegari & Rodriguez (2012), the equilibrium shape attained by the planet due to the gravita- 
tional interaction with the star, results in a prolateness, or equatorial ellipticity, which can be computed as (B — A)/C ~ 
(15/4)(M/m)(R/r) 3 . Knowing that C22 = (£/4)(B-A)/C we can obtain the magnitude of C22 for the planets with estimated 
radius. For the planets Kepler-lOb and 55 Cnce (see Section [3.2. II and 13. 2. 3p . we obtain C22 — 10 -3 and C22 — 2.4 x 10 -3 , 
respectively, assuming r = a currcnt . For sake of comparison with some rocky Solar System tidally evolved bodies, we note that 
C22 = 8.1 x 10" 6 in the case of Mercury (Smith et al. 2012 ), C22 = l.Ox 10" 5 for Titan (less et al. 2010) and C22 = 5.6 x 10" 4 
for Io (Schubert et al. 2004). 



3.2 Results 

The integration of the exact equations of motion usually requires a long computation time. For this reason, we decided to 
work with a small set of initial spin-orbit configurations and physical parameters. A statistical analysis taking into account 
many different initial conditions could give us a more accurate global picture of the past spin-orbit evolution of the planets. 
As pointed out by Correia & Laskar (2009), such statistical study is more adequate in the context of the averaged equations 
of the spin-orbit problem (see Section [4]). 



3.2.1 Kepler-lOb 

Kepler-lOb is a very close-in planet with a period of ~ 0.84 days, orbiting an old (~ 12 Gyr) Sun-like star (Batalha et al. 
2011). The physical parameters and the current orbital elements of the planet are summarized in Table 1. The estimated 
mass and radius values (m = 4.56m®, R = 1.4167?®) suggest a rocky composition of the planet, which can be considered as 
a super-Earth (Valencia 2011). The transit technique measurements do not allow a reliable estimation of the current orbital 
eccentricity, which was arbitrarily fixed at zero in this work. This choice is supported by our numerical simulations of the 
whole system, which includes also the exterior companion, Kepler-10 c. The results have shown that the long-term evolution, 
due to the coupled effect of tides and mutual interaction, could result in the complete circularization of both planetary orbits 
(Rodriguez et al., in preparation). 

We adopt the initial semi-major axis at aim = 0.02 au, which provides ei n i = 0.397 through equation ([7]). In addition, we 
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Figure 3. Time variation of the ratio Q/n for Kepler-lOb depending on the value of C%2- Different evolutions are observed, including 
captures in 2/1, 3/2 and 1/1 spin-orbit resonances (SORs). The solid black curve corresponds to the evolution for C22 = (i.e., when 
only the tidal torque exist) and coincides with the stationary tidal solution (STS, see text for discussion) of the rotation. The captures 
in the synchronous motion are detailed on the small panel. 



set f2ini = 2.7riini and £ = 0.35. For the tidal parameter, we assume Q' — 100 (a value usually adopted for rocky planets) and, 
solving equation we obtain feAt ~ 1.1 min. 

In Section [3.11 we have estimated that, for Kepler-lOb on its current orbit, C22 — 10 -3 . It should be emphasized that 
C22 would be smaller if the frozen- in shape of the planet was attained at larger distances than a currcnt . Hence, in order to 
investigate different possible scenarios, we consider three different values of C22, namely, 10 -3 , 10 -4 and 10 -5 . 

Figure [2] shows the variation of Q/n as a function of the eccentricity of Kepler-lOb. The curve with segments line 
corresponds to the stationary tidal solution (hereafter STS) of the rotation (0 = 0), obtained when only the tidal torque is 
considered. The STS is a function which depends only on the planet eccentricity, and is given by 

firtat/n = 1 + 6e 2 + 3e 4 /8 + 223e 6 /8 + 0(e 7 ) (9) 

(see equations (42) or (43) of Hut (1981); see also Section[4]for more details). The stationary rotation evolves with the tidally 
damped eccentricity and the system reaches synchronization (Q = n) only for circular orbits. Note that for small values of e, 
there follows that Q ~ n; for this reason, the STS is also called pseudo-synchronous solution (but this is valid only for e<l). 

The solid black curve in Figure [2] corresponds to the result of the numerical simulation performed with C22 = (only 
tidal forces are considered). In this case, the stationary rotation is attained almost instantaneously (when compared to the 
time of circularization). Note that both curves are in good agreement, at least for e < 0.4; hence, in this eccentricity interval, 
we can refer to the solution obtained with C22 = as STS. We also note that the curve for C22 = reaches the STS for Q/n 
close to two; however, this is incidental and unrelated with the 2/1 trapping. 

Figure. [3] shows the time variation of Q/n. For the three considered values of C22 > 0, the rotation rapidly encounters 
the value Q/n = 2 and thus the planet rotation is captured in the 2/1 spin-orbit resonance (hereafter SOR). On one hand, 
for C22 = 10 -3 and C22 = 10 -4 , the rotation escapes from the 2/1 trapping approximately at 80 Kyr, when the rotation 
drops to the value Q/n = 1, implying capture into 1/1 SOR and synchronous motion. However, it is worth noting that, for 
C22 = 10 -3 , the synchronism is attained in a form of jump, whereas for C22 = 10 -4 the rotation evolves without capture for 
a very short time before the 1/1 SOR (see small panel in Fig. [3|. We have also plotted the result for C22 = 0, showing that 
synchronization is attained around 400 Kyr. It is interesting to note that the introduction of the equatorial asymmetry of the 
planet reduces the time to reach synchronization; for instance, for C22 = 10~ 3 , it is reduced by a factor of 5. 

On the other hand, when C22 = 10~ 5 , the 2/1 SOR is destabilized at around 48 Kyr giving rise to the trapping in the 
3/2 SOR (Q/n — 3/2). Thus, the permanence of the rotation in the 2/1 SOR is temporary and its duration clearly depends 
on the C22 value. The rotation evolves in a Mercury-like state for some time, when a new instability occurs, after which an 
evolution free of trapping starts. Note that during the period of non-capture, the rotation returns to the stationary value, as 
can be seen in the dotted curve of Fig. [2] when e ^ 0.11. Finally, the synchronization is reached at around 270 Kyr. 

3 It is worth noting that the STS does not depend on any physical parameter of the interacting bodies. However, alternative tidal theories 
show that the STS would depend on the internal viscosity (see Ferraz-Mcllo 2012). 
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Figure 4. Time variation of semi-major axis (left panel) and eccentricity (right panel) of Kepler-10 b along the first 150 Kyr of evolution. 
The variation rates of a and e depend on the specific capture in which the rotation is trapped in. As different values of C22 lead to 
several captures, different scenarios are expected for the orbital evolution, despite the final values of the elements being the same for all 
C22 (see next figure). 
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Figure 5. Continuation of Figure [4] shows the final evolution. The orbital decay continues until circularization is achieved. The current 
semi-major axis of Kepler-10 b coincides with the one obtained in the simulation. 



Figures [4] [5] show the time evolution of semi-major axis and eccentricity of Kepler-lOb. We note that, when the rotation 
evolves under capture in resonance (that is, for C22 > 0), the rates of variations of the elements are larger than the one 
corresponding to the pure tidal case (C22 = 0), with the difference increasing as the eccentricity is damped (when the 
comparison is done for a given SOR). Moreover, for the same value of C22, the rate of orbital decay and circularization 
seems to depend on the specific SOR in which the planet rotation is trapped in. For example, in the case of C22 = 1CT 5 , the 
evolutions of a and e show two different regimes as the rotation evolves first in the 2/1 and then in 3/2 trappings (see Section 

EJ. 

We note that the time necessary to get an orbit with e = 0.05 is much greater for C22 < 10~ 3 . In fact, when the rotation 
evolves without capture, the time-scale for circularization increases, according to Figures [4] and [5] 

It is worth noting that, in the case of C22 > 0, the synchronous motion can be attained even for e / 0. For instance, 
in the case of C22 = 10~ 4 , the system is captured in 1/1 SOR at roughly 120 Kyr, when e ~ 0.079. This fact confirms the 
classical result that, for eccentric orbits, the synchronous motion is only possible if there exists an additional torque which 
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Figure 6. Evolution of the orbit of Kepler-lOb in the plane (a, e). The solution of equation (0 is shown by the dashed line and is 
compared to the result of the numerical simulation for the case C22 = 10 — 3 , giving a good agreement (the curves for the other two C22 
values also agree with the solution of equation J7J, but are not shown here for sake of simplicity). The starting position of the planet is 
indicated by the star symbol. 



counteracts the tidal torque (see Goldreich & Peale 1966; Ferraz-Mello et al. 2008 and references therein). We also note that 
after circularization, the 1/1 SOR is the only possible resonant motiorQ. 

Figure [6] shows the evolution of the system in the plane of orbital elements a and e. The solid line shows the result of 
numerical simulation in the case of C22 = 10 -3 . The dashed line corresponds to equation flJJ, which is independent on the 
value of C22 - Note as both curves are in very good agreement, in such a way confirming our assumption about the conservation 
of the orbital angular momentum (Section 13. II) . 

3.2.2 G J 3634 b 

GJ 3634 b is a super-Earth planet with mass m = 7.0m® orbiting a M2.5 dwarf each 2.65 days with an eccentricity of e = 0.08 
(Bonfils et al. 2011). Since there is no transit detection for this planet, its radius is unknown. However, Valencia, Sasselov & 
O'Connell (2007) have shown that there is a maximum radius that a rocky terrestrial planet may achieve for a given mass. 
A larger radius would imply a large amount of water and even the possibility of a massive volatile envelope. Indeed, for a 
planet with mass 7.5m® the maximum terrestrial radius would be 11,600 km (~ 1.827?®), which result in a mean density of 
p ~ 1.25p®. Hence, we assume for GJ 3634 b that p — p® to obtain R ~ 1.97?®. 

The adopted initial orbital elements are ai nl — 0.033 au and e- m i = 0.369, while fiini = 3.7ni n i, £ = 0.35 and Q' = 100 
(feAt = 2.1 min). At the current value of the semi-major axis the estimated value of C22 is C22 — 1.6 x 10 -4 (Section 13. 1[) ; 
in this work, we perform the analysis using two values of C22: 10 -4 and 10 -5 . 

Figure [7] shows the variation of 0,/n as a function of e (left panel) and the time variation of a (right panel). For C22 = 10 -4 , 
the rotation is rapidly trapped in the 2/1 SOR, remaining in such resonant configuration during the most time of the tidal 
eccentricity damping. Only at low eccentricities, close to 0.032, the system leaves the 2/1 SOR dropping directly into the 
synchronous motion. 

At variance, for C22 = 10 -5 , the rotation quickly reaches the stationary (non-resonant) state given by the STS curve 
(C22 = 0) in Figure [7] (left). The system evolves along the STS curve up to close approach to the 3/2 SOR. Then, the system 
is captured and remains in this resonance until dropping in the synchronized state, at e ~ 0.012. 

In both cases, the final configuration of the planet's orbit is characterized by circularization and synchronization, as in 
the case of Kepler-10 b. However, the current eccentricity of the GJ 3634 b is different from zero (its location at 0.08 is shown 
by the dotted vertical line in the left panel of Figure [7} . Note that the current eccentricity is reached when the rotation still 
evolves in the 2/1 or 3/2 SOR (depending on C22). This fact suggests that GJ 3634 b could be presently involved in one of 
such resonant trappings, in particular, it could be found in a Mercury-like rotational state (Q/n ~ 3/2). However, it should be 
kept in mind, that this conclusion depends on the choice of specific values of C22 and Q' . Indeed, the results of an additional 
simulation performed with C22 = 10 -6 , have shown that the trapping of the system in the 3/2 SOR is destabilized at e ~ 0.11. 
Since e CU rrcnt = 0.08 < 0.11, the planet rotation would be currently evolving without capture (along the STS) for such a small 
C22 value. 

The semi-major axis evolution (Figure [7] right panel) shows that the tidal decay is faster when the rotation of the planet 
is trapped in the 2/1 and 3/2 spin-orbit resonant motions. The fact that the current eccentricity of the planet orbit is different 



4 Except if an additional mechanism, such as impact events, acts to destabilize the synchronous motion (see Wieczorek et al. 2012; 
Correia & Laskar 2012) 
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Figure 7. The variation of Q/n as a function of e (left panel) and the time evolution of a (right panel) for the planet GJ 3634b. 
Depending on C22, the planet rotation may be currently evolving in either 2/1 or 3/2 SOR. Note that the time-scale of orbital decay 
is drastically reduced provided the planet can support a permanent deformation resulting in a nonzero C22 value. The vertical and 
horizontal dotted lines indicates the current values of e and a, respectively 



from zero suggests that, in the absence of other perturbations, the orbital decaj(f] of the planet is presently in progress and, 
consequently, the final semi-major axis will be smaller than the current one. 



3.2.3 55 Cnce 

55 Cnce is a hot super-Earth planet, with the mass m = 8.63m ffi and the radius R = 2.17i?©, orbiting a Sun-like star each 
0.74 days with e = 0.057 (Winn et al. 2011; Gillon et al. 2012). The averaged density of the planet is comparable to that of 
Earth, suggesting a rock-iron composition supplemented by a significant mass of water, gas, or other light elemental (Winn 
et al. 2011). 

For numerical simulations, we adopt a-mi = 0.0186 au, ei n i = 0.4, fii n i = 2.7rti n i and £ = 0.35. The tidal history of the 
system indicates that the magnitude of Q' must be higher than the commonly adopted values for rocky planets (Ferraz-Mello 
2012). In this work, we use Q' — 5.5 x 10 3 ; this value allows us to obtain A^At ~ 1.1 s. We perform numerical simulations for 
two different values of C22: 10~ 3 and 10~ 4 . For the purpose of further comparison with the results reported in Callegari & 
Rodriguez (2012), we also consider the values C22 = 2.45 x 10~ 3 and C22 = 2.45 x 10 -4 . 

Figure [S] shows the evolution of rotation as a function of the planet's eccentricity. The captures inside the 5/2, 2/1, 3/2, 
4/3 and 1/1 resonances can be observed. For C22 = 2.45 x 10 -3 , the planet rotation evolves under capture in the 3/2 SOR 
almost up to the orbital circularization, with fl/n roughly varying between 1.46 and 1.54. 

In a recent work, Callegari & Rodriguez (2012) have shown that, for e = 0.057 and (B ~ A)/C = 0.028 (or C22 = 
2.45 x 10~ 3 ), the phase space of the rotation (0/n,0) of 55 Cnce, between the main spin-orbit resonances, is dominated by a 
large region of chaotic motion. This region increases as larger values of the eccentricity are considered. 

Capture in spin-orbit resonance can work as a mechanism to prevent a chaotic rotation, provided that the dissipation is 
large enough in order to ensure the corresponding trapping. According to our simulations (see Figure [S]), destabilizations of 
different resonances occur as the eccentricity is tidally damped (see also next section), giving rise to new temporary trappings 
with a subsequent evolution to synchronism. In this way, a chaotic behavior of the rotation can be avoided by a continuous 
evolution under resonant motion. It is important to note that, at least for the cases shown in Figure [5] there is no interval of 
non-resonant motion between two successive trappings, as in the case of Kepler-lOb (see Figure[3]for C22 = 10~ 4 and 10 -5 ). 

An additional simulation done with C22 = 2.45 x 10 (not shown in figures), has resulted in a 2/1 trapping until e ~ 0.04, 
indicating that the planet rotation may be currently evolving in resonant motion, for such a small value of C22- 

Figure [5] shows the time variation of semi-major axis. It is worth noting that the orbital decay results in a more rapid 
evolution towards o OU rrent when the planet rotation is trapped in a SOR. Also, the rate of variation depends on the specific 
resonance. From Figure we see that the evolution when C22 ~ 10~ 4 , (capture in 2/1) is almost two times faster than the 
one for C22 = 2.45 x 10~ 3 (3/2 capture). We will return to this discussion in the next section. 



5 In fact, a oc e 2 when fl = n or Q ~ n (see equation H20l l). 

6 The physical composition and the suggested value of Q' (5.5 
instead of a super-Earth should not be ruled out. 



10 3 ) indicate that the possibility of 55 Cnce to be a mini-Neptune 
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Figure 8. The variation Cl/n as a function of e (in logarithmic scale), for the planet 55 Cnce. We have used C22 = 10 in addition to 
the value adopted in Callegari & Rodriguez (2012) C22 = 2.45 X 10 — 3 . Also values ten times smaller were taken into account for sake 
of completeness. The vertical dotted line indicates the current value of e, namely, 0.057. The trappings in the RSOs 5/2, 2/1, 3/2, 4/3 
and 1/1 appear during the evolutions. Note the large amplitude oscillation of f!/n inside the 3/2 and 4/3 cases with C22 of the order of 
10 -3 . Depending on C22, 55 Cnce can be an example of a planet currently evolving in the 5/2, 2/1 or 3/2 SOR. 
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Figure 9. Time variation of semi-major axis for 55 Cnce. The curves show the evolution for four values of C22 together with the STS 
curve. As for previously analysed planets, the orbital decay is faster when the planet rotation is trapped in a resonant motion. The 
current value of a is obtained for all values of C22 ■ 



4 AVERAGED EQUATIONS 

In order to better understand the results previously obtained through numerical simulations of the exact equations of motion, 
in this section, we develop an analytical approach in the context of the averaged equations of the spin-orbit problem. 

The equation of motion for the spin in equation © cannot be solved analytically. However, because we are interested in 
the evolution inside spin-orbit resonances, we will consider planetary spin rates having values such that Q ~ pn, where p is 
a half-integer. Our investigation is restricted to the cases p = 1, p = 3/2 and p = 2 but can be easily extended to any other 
SOR. 

We introduce the "resonant" angle S — 8 — pi, where I is the mean anomaly of the planet orbit and 9 was defined in 
Section [2] Note that the angle S gives the orientation of the long axis whenever the planet passes through the pericenter of 
its orbit (i.e., 5 — 8 when I = 0, see Figure [TJ. It is easy to show that, inside the spin-orbit resonance given by the parameter 
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p, the resonant angle is a slowly varying variable when compared to the orbital period of the planet. Indeed, we can write 
that 9 = 5 + pn, where 6 — Q, ~ pn and, consequently, 5 ~ 0. Now, replacing 8 by 5 + pi in equation ((3J and averaging over 
the orbital period, we obtain the classical result for the averaged torque on the planet produced by the permanent equatorial 
deformation (see Goldreich 1966; Goldreich & Peale 1966): 

< T 22 >= -6mR 2 C 22 n 2 H(p, e) sin 25, (10) 

where H (p, e) are power series of e with the leading term of order e 2 ' p_1 ' (see Appendix A for a literal expansion up to seventh 
order) . 

The averaged tidal torque on the planet is given by 

< T tide >= if tide Q T(p,e), (11) 

where Ktide = 3k 2 AtMR 5 and T{p,e) is a suitable function, whose expansion in powers of e is shown in Appendix A (Hut 
1981). 

At resonance, the equilibrium solution of the rotation is given by the condition 
Ctl =< T tidc > + < T 22 >= 0, (12) 
and, resolving for C 22 sin 25, we obtain 

C7 22 sin2^I fc2 At^(^ 3 ^4. (13) 
2 m \ a J H{p, e) 

If we know the maximum value that 5 can reach without disturbing the trapping in the p-SOR, we may obtain the minimum 
value of C22 necessary for capture in the same resonant motion. 

4.1 1/1 Spin-Orbit Resonance 

Replacing p = 1 in equation (|13l) and expanding the ratio g^ej in powers of e, up to order 6, we obtain 

C 22 sm25 = — k 2 Atn— I - (48e 2 + 483e 4 + 2624e 6 ) + C(e 8 ). (14) 
16 m \ a J 

We note that, for 0(e 4 ), the above expression is the same obtained in Ferraz-Mello et al. (2008), for synchronous asymmetric 
planets 0. 

4-1.1 Semi-major axis evolution 

The averaged variation of the orbital elements can be obtained using the Euler-Gauss equations (see Beutler 2005, p. 230), 
which are written in terms of the radial and transverse components of the perturbing force. The result for the semi-major 
axis is 

< d 22 >= - 3nfi2 ^ Sin ^ (-288 + 720e 2 - 234e 4 + 35e 6 ) + 0(e 8 ). (15) 
Replacing 6*22 sin 2(5 by the corresponding expression given in equation (|14[) . we obtain 

< a 22 >= 3fc 2 At^ (|) 5 n 2 a (l2e 2 + ^e 4 + ^i e 6 ) + 0(e 8 ). (16) 

In the context of the linear tidal model, the averaged variation of the semi-major axis due to tides is given by equation 
(9) of Hut (1981), valid for any values of p and e. The expansion up to C(e 8 ) reads 



+ 0(e 8 ). (17) 



< atidc >= 3k 2 At— ( — ) na 2(-l+p) + (-46 + 27p)e 2 + -(-480 + 191p)e 4 + -(-13530 + 3961p)e 6 

Til \ a J L 4 8 

Replacing p = 1 in the above equation we obtain 

■xt. a + M f R \ 5 2 Ciq 2 , 867 4, 9569 «A , 8 . , - 

< atido >p=i= -3k 2 At— I —J n a I 19e + — — e H — e J + U(e ). (18) 

The total averaged variation of the semi-major axis is then given as 

< a > p =i = < dtide >p=i + < a 22 > (19) 
and thus 

7 In this case, the association e 2 = (2f! — n)At = nAt should be done, where t 2 is the phase lag of the tidal wave whose frequency is 
(2fi - n). 
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< a > p=1 = -3k 2 At™ (*) n*a he 2 + ^ + + G(e»). (20) 

m \a J \ 4 4 J 

At the same order of approximation, the above equation coincides with the result of previous works (see Ferraz-Mello et al. 
(2008) and references therein). 

When only the tidal torque on the planet is taken into account, the non-resonant STS is given by equation ([9]). In this 
case, the synchronism is possible only for circular orbits. If instead of p — 1, we replace p = p s t a t = f2stat/n in equation (|17[) . 
we obtain 

< a tidc > P=Pstat = -3foAt^ ) 5 n 2 a he 2 + 54e 4 + ^e 6 ) + 0(e s ). (21) 

We define r„ as the ratio between averaged semi-major axis variation in the p-SOR (p = 1, 3/2, 2, ...) and in the tidal 
stationary solution of the rotation (p = p s tat), that is 

' " : " . (22) 



< fitidc >p=p s tat 

Comparing equations (|20[) and (I21|l . we note that r^ 1 > 1, indicating that the rate of semi-major axis decay is stronger when 
the rotation is locked in the 1/1 SOR, compared to the case when it evolves following the STS. This result is in agreement 
with Wisdom (2008), who has shown that the rate of energy dissipation in a body with synchronous rotation is larger than 
in a asymptotic non-synchronous rotation. 

4-1.2 Eccentricity 

We repeat the procedure described in the previous section to obtain the averaged variation of the planet eccentricity. The 
resulting expressions are: 

< e 22 >= -9fc 2 At^ (^) 5 n 2 (e 3 + i^e 5 ) + Q(e 7 ), (23) 

< e tidc > P =i= -h 2 At^- (*V n 2 he + 113e 3 + ^e 5 ) + 0(e 7 ). (24) 

< etidc > P = Pstat = --kiAt—l-) n 2 (28e + 188e 3 + 917e 5 ) + C(e 7 ). (25) 

< e > p= i= --k 2 At— (- n 2 (28e + 476e 3 + 2825e 5 ) + C(e 7 ), (26) 

Similarly to the semi-major axis variation case, rf =1 > 1, where r£ is the analogous for the eccentricity of equation (1221) . 



4.2 3/2 and 2/1 Spin-Orbit Resonances 

We analyse in this section the spin-orbit behavior when the planet evolves in the 3/2 and 2/1 SORs, employing the same 
method described in the previous sectiorjf]. We have 



/•^ o, A ,M fR\ 5 2 /l 5 2 774 4 8019 6 \ 8 . 

< a > p=3/2 = Sk 2 At- [- J n a - -e + — e + — e j + 0(e ), (27) 
and 



< a > P=2 = ~3fc 2 At^ (|) n 2 a (2 - 2e 2 - 3e 4 + ^e 6 ) + 0(e 8 ). (28) 

Figure [TOl shows a plot of for the three analysed resonances as a function of the eccentricity. On one hand, we note that 
r^ 2 increase for small eccentricity, indicating that the major effect on the orbital evolution due to the resonance 
captures occurs as e decreases, in agreement with the results of the numerical simulations (see Figures 2] [7] and . On the 
other hand, for large e, r a is small (mainly in the case p = 3/2) as also shown in Figures [4] and [7] where there is a small 
difference between the evolutions with and without C22 (i.e, r a — 1). 

When the synchronous motion is analysed, we note in Figure [TD] that rg =1 is an increasing function of e, indicating that 
the orbital decay when the rotation is trapped in the 1/1 SOR is large for large e. However, as simulations show, synchronism 
is attained for very low values of e, making the orbital decay almost negligible. 

The results above are also valid for r\. It is interesting to note that r v a and rf are only functions of the eccentricity and 
do not depend on any physical parameter. 

It is worthwhile mentioning that T 22 averages to zero unless the condition ~ pn is assumed (whit p integer), that is, 
the averaged orbital evolution due to the triaxial torque is expected only for resonant cases. 

8 We stress that H (p > 1, e) has no term independent of e and thus equation II13I I cannot be expanded around e ~ 0. 
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Figure 10. The variation of rj as a function of orbital eccentricity, for p = 2, 3/2 and 1. In the cases of p > 1 the ratio of rates of 
variation of a in the SOR and in the STS increases for small e, showing that the orbital decay is enhanced (when compared with the 
evolution in the STS) whenever the planet rotation evolves under resonance with the orbital motion. 



4.3 Critical C22 

As it has been noted, equation (|13[) can be used to obtain the minimum value of C22 necessary for the stability of the p-SOR. 
Calling C22 the one value given by the right hand of equation (|13[) . the stability condition for a resonant spin state requires 
that C22 > C2V The above condition is frequently found in literature as strength criterion of the resonance (see Murray & 
Dermott 1999, p. 205). 

Figure fTTI shows the functional dependence of the absolute value of C22 with the eccentricity, obtained for Kepler-lOb. 
The dashed horizontal lines represent three values of C22 used in the numerical simulations shown in Section 13.2.11 namely, 
10" 3 ,10" 4 and 10~ 5 . 

The case p = 2, shown by the blue curve, indicates that C22 = 10~ 5 for e ~ 0.3, C22 = 10~ 4 for e ~ 0.1 and CJ2 = 10 -3 
for e ~ 0.03. A good agreement is obtained doing the comparison with the results of the numerical simulations for Kepler-lOb, 
shown in Figure [T2l Note as the 2/1 SOR destabilizes when e ~ 0.28, e ~ 0.1 and e ~ 0.03. 

The same can be observed for the 3/2 SOR, in which the rotation escapes from the trapping when e ~ 0.11 if C22 = 10 , 
in concordanc^EI with the result of Figure [TT] (orange curve). 

Since for roughly e ^ 0.1 the condition C22 > C22 is satisfied more easily for p = 1 as compared to p > 1, there is no 
problem in understanding the occurrence of synchronous motion for very small eccentricity, as can be also verified in Figure 

cm 

It is worth noting that instead to compute a critical C22 value, the condition for stationary solutions < Ttido > + 
< T22 >= CQ = also enables us to obtain a critical value of e for a given fixed C*22- The resonant motion would destabilize 
whenever the eccentricity becomes smaller than the critical value for each SOR. This approach is more convenient in those 
cases for which the observations allow an estimation of C22 (see Correia & Laskar (2004) for the study of the rotation of 
Mercury). 



5 DISCUSSION AND CONCLUSIONS 

The numerical exploration of the planet rotation which is subject to the conservative triaxial gravitational torque of the 
star shows different behaviors, including the oscillation around the equilibrium points of spin-orbit resonances (see Callegari 

6 Rodriguez 2012 and references therein). When dissipative effects produced by tidal torques are included, the rotation of 
close-in bodies can be captured in a resonant motion (Goldreich 1966; Goldreich & Peale 1966; Celletti & Chierchia 2008). 
The capture in a specific resonance depends on the eccentricity e and on parameters related to tides and equatorial ellipticity 
(Q and C22, respectively). Hence, since e is tidally damped during the inward migration, each specific capture should become 
unstable and the rotation can achieve another resonant state, which, at the same time, should result in a temporary trapping. 
When the orbit is circularized, the final evolution results in the synchronization between the orbital and rotational periods 
(i.e, the 1/1 spin-orbit resonance). 

In this paper, we analysed the coupled spin-orbit evolution of Kepler-lOb, GJ 3634b and 55 Cnce short-period super- 
Earth planets. The numerical simulation of the exact equations of motion have shown that, depending on C22, the rotation 



9 Small differences between the result of numerical simulations and the analysis of averaged equations may arise because for the con- 
struction of Figure [TT] we use a = a currcn t, taking also into account that truncated series in e were used for the same figure. 



14 A. Rodriguez et al. 




Figure 11. Critical value of C22 given by equation J 1 3 1) and calculated with parameters of the planet Kepler -10 b. The calculation of 
can be used as a guide to investigate the stability of a given SOR as a function of eccentricity (see text for discussion). Note that 
in the case of p = 1 we have C| 2 = for e = (see equation 1 1141) ) , indicating that the synchronism is the only possible resonant motion 
for circular orbits (see Goldreich & Peale 1966; Celletti & Chierchia 2008; Callegari &; Rodriguez 2012). 




0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 



e 

Figure 12. The variation of Q/n as a function of the orbital eccentricity for Kepler-lOb, resulting from the analysis of numerical 
simulations of the exact equations of motion (see Section |3. 2D . If a comparison with the results of Figure [TT] is done, we obtain a good 
agreement for the values of e for which each SOR destabilizes (see text for details). 



displays a number of resonance trappings for which Q ~ pn, including the cases p = 5/2, 2, 3/2, 4/3 and 1. Moreover, 
the results indicate that GJ 3634 b can be an example of an exoplanet with a rotation state similar to Mercury's rotation 
(p — 3/2), provided the body figure can support a permanent equatorial asymmetry as small as C22 = 10 -5 . 55 Cnce can 
currently evolve under 5/2 or 3/2 trappings if C22 ~ 10 -3 and under 2/1 if C22 ~ 10~ 4 or even if C22 is as small as 2.45 x 10~ 5 . 
The rotation of 55 Cnce may have a chaotic behavior for C22 ~ 10~ 3 (Callegari & Rodriguez 2012). However, we have shown 
that a scenario is possible in which the rotation evolves to synchronism through successive trappings in SOR. Indeed, the 
eccentricity damping due to tides provides the destabilization mechanism required to escape from one given capture and 
subsequently evolve to lower order resonances (see Correia & Laskar 2004, 2009). 

We have shown that different spin-states corresponding to various resonance captures can occur in extrasolar systems. 
Although the 1/1 spin-orbit coupling is the final state of evolution, it is not reached in every case as present configuration. For 
current non-circular orbits, trappings in the 3/2 SOR or other resonance configurations can be sustained over long times. The 
latter would depend on the orbital eccentricity, the planet's C22 value and its dissipation factor. Given that these parameters 
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are not well-determined for terrestrial extrasolar planets, we note that the 1/1 coupling may not be the current rotation state 
of some planets. This would have consequences on the stellar insulation and therefore on the evolution of the atmosphere and 
surface of the planet. 

When only the tidal torque is taken into account (C22 = 0), the rotation evolves to the stationary tidal solution (STS). 
We have shown that under assumption of C22 > 0, the evolution along spin-orbit resonance (SOR) leads to higher rates of 
orbital decay and circularization than the corresponding ones when the rotation evolves in the STS. Moreover, for a given C22 
value, the variation rates of the elements depend on the specific trapping in which the rotation is trapped in, and increases 
with p (see Figure [3J. 

We have also calculated the averaged rates of variations of semi-major axis and eccentricity including the gravitational 
and tidal torques. Moreover, through the analytical calculation of the critical value of C22 it is possible to predict the values 
of the eccentricity for which a given SOR would become unstable. The comparison with the direct numerical simulations have 
shown a good concordance for the analysed planets. 

The importance of the C22 value on the rotational history of the planet is illustrated in the example of Kepler-lOb. It 
is worth noting that the planet rotation of Kepler-lOb can evolve under the influence of the 2/1 SOR for longer time-scales 
(i.e., small e) provided the planet figure can support a permanent equatorial deformation characterized by C22 ~ 10 -3 . For 
smaller C22, the 2/1 resonant state would become unstable for relatively large eccentricity (Figures 1111 and [T2l . 

We stress that the results of the numerical simulations provide different possible spin-orbit evolutions depending on the 
considered initial values of elements and rotation. The Q value and the specific tidal model would also play an important role 
in determining the rotation history of the planet, since the temporary trapping in SOR depends on the rate of eccentricity 
damping, which is Q-dependent. In fact, in a recent paper, Makarov (2012) has applied the model by Efroimsky (2012) and 
successfully explained the current resonant state of Mercury within the considered model. 

Our results concerning the association of the critical C22 value with the capture probability are not new in the context 
of averaged equations of the spin-orbit problem (see for instance Correia & Laskar 2004, 2009). However, in this work, we 
are able to compare the results of the averaged equations with those resulting from the numerical simulations of the exact 
equations of motion. 

A model of planet interior including the effect of viscous friction at the core-mantle boundary would provide quantitatively 
different results (see Peale 2005, and Correia & Laskar 2009, for the case of Mercury). According to previous works the chances 
of resonant capture increase if the latter effect is taken into account (Goldreich & Peale 1967; Peale & Boss 1977). However, 
in order to highlight the contribution of C22 to the orbital evolution, in this paper we adopt a simple model neglecting a 
possible layered structure of the planet. 

In the present study we have assumed constant C22 values. However, depending on the adjustment to hydrostatic equi- 
librium C22 may be itself time-dependent. The time-scale for adjustment to the completely relaxed state depends on material 
parameters, e.g., rigidity and viscosity which are highly temperature dependent. Furthermore, it will depend on the spin-orbit 
coupling the planet is trapped in. For the application to terrestrial planets given here the constant, i.e. frozen, C22 configura- 
tion is justified because adjustment of solid rocky planets to equilibrium usually occurs on geological time-scales. Therefore, 
the states described here do not necessarily represent hydrostatic configurations. 

In this work we have considered a model of one interacting planet, despite the fact that Kepler-lOb and 55 Cnce belongs 
to multiple-planet systems with exterior (more massive) companions. Previous works have shown that the mutual secular 
interaction with the companion may affect the eccentricity evolution of the tidally evolved planet, increasing the circularization 
time-scale (Mardling 2007; Rodriguez et al. 2011a; Rodriguez et al. 2011b). Moreover, the planetary perturbations can produce 
large excursions of the innermost planet eccentricity, providing a mechanism for increasing or decreasing the capture probability 
of the spin-orbit resonances (Correia & Laskar 2004, 2009). However, a numerical simulation of 55 Cnc system including the 
exterior companion 55 Cncb has shown little differences with the single planet case. Indeed, the evolution of the inner orbit 
eccentricity is almost not affected by the mutual interaction between the planets and thus the rotation evolution of 55 Cnce 
can be safely analysed within the one planet model. An investigation about the rotation of asymmetric close-in super-Earths 
including an exterior companion will be addressed in a future investigation. 
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APPENDIX A 

We show here the explicit expressions for the functions H(p, e) and T(p, e). The exact form of H(p, e) is given by 




(29) 



The expansion up to order six in e is displayed in Table 2 for p = 1,3/2 and 2. 
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Table 2. Expansion of the function H(p,e) given by equation J29H for 0(e 7 ) for the principal spin-orbit resonances. 
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1 6 2 , 13 _4 _ 35 6 
2 ' 16 288 



3/2 



1 P _ 123^3 , 489 5 
2 16 ' 128 



2 



12 P 2 _ 115 ,4 , 601-1 
2 6 T 48 



The exact expression for the function T(p,e) is given by (Hut 1981) 



T(p,e) = (1 - e 2 r° [/2(e) ~ (1 - e 2 ) 3/2 /5(e)p] , 



(30) 



where 



He) 




5 , 
— e 
16 



(31) 



and 



Me) 



i i o 2 . 3 4 

1 + 3e + -e . 



(32) 
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